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Abstract -We analyze nonlocally coupled networks of identical chaotic oscillators with either 
time-discrete or time-continuous dynamics (Henon map, Lozi map, Lorenz system). We hypoth¬ 
esize that chimera states, in which spatial domains of coherent (synchronous) and incoherent 
(desynchronized) dynamics coexist, can be obtained only in networks of nonhyperbolic chaotic 
systems and cannot be found in networks of hyperbolic systems. This hypothesis is supported by 
analytical results and numerical simulations for hyperbolic and nonhyperbolic cases. 


Introduction. — Chimera states in dynamical net- 
i__iworks of nonlocally coupled chaotic oscillators have re¬ 
cently attracted much attention [^. They represent an 
' ^ intriguing phenomenon where an ensemble of identical 

_elements with symmetric coupling spontaneously splits 

into spatially separated coexisting domains of coherent 
00 (synchronized) and incoherent (desynchronized) dynam- 
in ics. Since their first discovery in systems of coupled phase 
^^oscillators they have been found in a broad range of 
diverse models , ranging from time-discrete maps 
and time-continuous chaotic models , Stuart-Landau os- 
dilators |7[|^, globally coupled lasers , Van der Pol os- 
. .dilators [l2[|13 FitzHugh-Nagumo neural systems pp^ , 
^ ^ population dynamical models 16 to autonomous Boolean 
networks 


17 


They have recently also been observed ex- 

chem- 

electronic 


5 _j perimentally, e.g., in optical light modulators 18 

epical 


19, mechanical ^ 21 


22 23 , optoelec¬ 


tronic 24 , and electrochemical [25|[26] oscillator systems. 
In real-world systems chimera states might play a role, 
e.g., in the unihemispheric sleep of birds and dolphins 
in neuronal bump states 28 29 , in epileptic seizure 


in power grids 31 , or in social systems 32 


27 


30 


In spite of intense theoretical study, no universal mech¬ 
anism for the formation of chimera states in different 
systems has yet been established. It remains unclear 
what the common features of the different specific mod¬ 
els with different local dynamics are. Chimera states 
emerge as a hybrid state on the transition between com¬ 


pletely synchronized coherent and completely desynchro¬ 
nized incoherent states with decreasing coupling strength, 
and such coherence-incoherence bifurcation scenarios have 
been found to display unversal features in many different 
time-discrete and time-continuous models iElIl- 


In this work we analyze networks of identical chaotic 
oscillators with nonlocal coupling, with the aim to under¬ 
stand the mechanism of the emergence of chimera states 
in different systems. We study networks of both time- 
discrete chaotic maps (Henon and Lozi map) and time- 
continuous three-variable chaotic systems (Lorenz model). 
In the present work we observe that the emergence of 
chimera states appears to be linked with the property 
of nonhyperbolicity of the dynamical system, while hy¬ 
perbolic system^ appear not to exhibit chimeras. In a 
number of papers chimeras have been reported in ring 
networks of chaotic nonhyperbolic systems exhibiting a 
period-doubling route to chaos, e.g., the logistic map [^, 
the cosine map 18 and the Rossler attractor [^, while 
chimeras of oscillating states have not been observed for 
the Lorenz attractor 33 , which is hyperbolic for the stan¬ 
dard parameters. In particular we study the Henon map 
34 , which is an example of a nonhyperbolic dynamical 
system, like the Rossler attractor, and the Lozi map as 
an example of a hyperbolic system. By tuning the pa¬ 
rameters of the Lorenz model we show that chimeras arise 


^ These systems are reputed hyperbolic because of non-existence 
of stable regimes and its bifurcations. 
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in the nonhyperbolic regime, while they are suppressed in 
the hyperbolic regime. 

Hyperbolic and nonhyperbolic chaotic systems. 

— Fist, as an example of a nonhyperbolic dynamical sys¬ 


tem, we study the Henon map 34 

= 1 - a{xy + y\ 




= 


( 1 ) 


where t is the dicrete time, x*, y* are the dynamical vari¬ 
ables, and a > 0 and /3 > 0 are parameters of the map. It 
should be noted that the case of 0 < /3 < 1 corresponds 
to a dissipative diffeomorphism. The Henon map under 
high-ratio compression (/3 —)• 0) reduces to the logistic 
map = 1 — a{x^)^. It belongs to the wide class of 
the Feigenbaum maps, characterized by a quadratic max¬ 
imum. When the parameters are changed, the map Q 
demonstrates a period-doubling bifurcation cascade and 
periodic orbits which are in agreement with Sharkovsky 
ordering [35) . 

The Henon map can be used for a qualitative descrip¬ 
tion of bifurcation phenomena in three-dimensional time- 
continuous systems of spiral-type chaos because of the 
saddle-focus separatrix loop. As shown, for example, in 


36 the spiral-type chaos in a three-dimensional system 
gives a two-dimensional map in the Poincare cross section, 
equivalent to Eq.Q. The Henon map Q is an example of 
a nonhyperbolic dynamical system with bifurcation of ho¬ 
moclinic connections of stable and unstable saddle point 
separatrices 37 . This bifurcation is the basis for the emer¬ 


gence of the nonhyperbolic attractor. It means that the 
Henon map is the simplest model of a chaotic map which 
can be obtained as the Poincare section of the nonhyper¬ 
bolic class of Rossler oscillators. Figure [^a) illustrates 
the effect of homoclinic tangency of the stable (IF®) and 
unstable (IF“) saddle point manifolds in the Henon map. 
The homoclinic tangency leads to the complex structure 
of the Henon attractor shown in Fig.[^b). The nonhyper¬ 
bolic attractor of the Henon system is characterized by the 
coexistence of several regular and chaotic attracting sets. 
Their attractor basins (white and purple regions) have a 
fractal structure. Upon change of parameters, these sets 
undergo various bifurcations, which are accompanied by 
drastic changes in the structure of the attracting set. 

The time-continuous Lorenz system belongs to another 
class of dynamical systems. It is an example of a hyper¬ 
bolic attractor without regular stable subsets like equilib¬ 
rium states or limit cycles. The Lorenz attractor is re¬ 
silient to slight changes in parameters and perturbations 
of the system equations. The Lozi map 38 belongs to 


the class of chaotic maps which can be obtained as the 
Poincare section of the hyperbolic class of Lorenz-type at¬ 
tractors: 


x^'^^ = 1 — a\x*\ + y*, y*^^ = Px*. (2) 

The stable and unstable saddle point manifolds of the Lozi 
attractor intersect transversally (without tangency). The 




X 
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Fig. 1: (a) - Stable {W”) and unstable (IF“) manifolds of a 
saddle point q of the Henon map for a = 1.3 and /3 = 0.3. 

(b) - Attractors in the Henon map and their basins of attrac¬ 
tion in the phase plane {x,y) for a = 1.078 and P = 0.3. 
The white region is the basin of attraction of the four-band 
chaotic attractor (black dots), the purple (dark grey) region 
corresponds to the basin of attraction of the six-band chaotic 
attractor, and trajectories from the yellow (light grey) regions 
diverge to infinity. 
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Fig. 2: (a) - Stable {W”) and unstable (IF“) manifolds of a 
saddle point q in the Lozi map for a = 1.7 and /3 = 0.3. 

(b) - Lozi attractor Go and the basin of its attraction Gi 
(white) for a = 1.5 and P = 0.3. Trajectories from the yel¬ 
low (light gray) region diverge to infinity. 


property of transversality is resilient to parameter varia¬ 
tions. This means that the hyperbolic Lozi attractor is 
robust (structurally stable). It is illustrated in Fig. [^a). 

Figure [^b) shows that the Lozi attractor is a unique 
connected attracting set with bounded and homogeneous 
attractor basin in the phase plane. The attractor topology 
is robust to slight changes in the control parameters. 

Dynamics of ring networks of nonlocally coupled 
Henon and Lozi maps. — We consider a ring of N 
nonlocally coupled two-dimensional maps, described by: 

i+P 

( 3 ) 

= Pxl, 

where t is the discrete time, N is the number of elements 
in the ring, i = 1,2... N is the index of the element, all 
indices are modulo N, a is the coupling strength, P is 
the number of neighbours on either side, and f{x,y) = 
1 — -b y (Henon map, Eq.Q), or f{x,y) = 1 — a\x \ + y 
(Lozi map, Eq.Q). We introduce the parameter r = P/N 
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as the coupling radius. 

First, let us analyze the nonlocally coupled Henon maps 
0. Figure 0 illustrates all main regimes of this sys¬ 
tem. It shows that the system of coupled Henon maps 
demonstrates the same regimes as coupled logistic maps 

. Their (r, a) parameter planes have the same quali¬ 
tative and quantitative structure with coherence tongues 
(H, B, C). The area E in Fig. 0 corresponds to the regime 
of the spatially incoherent states, D is the region of the 
completely synchronized chaotic states. The regions 
B, and C correspond to the spatially coherent states with 
wave numbers fc = 3, fc = 2, and k = 1, respectively. 



Fig. 3: Main regimes of nonlocally coupled Henon maps (a) 
and Lozi maps (b) in the (r, cr) parameter plane. (a) - 
Regimes of coherence with wave numbers k = 1,2,3 (regions 
C, B, A, respectively). The horizontal line at Uc ~ 0.5 
marks the coherence-incoherence transition, (b) - White re¬ 
gion corresponds to multistable regime of stationary patterns 
and travelling wave profiles (insets show profiles with differ¬ 
ent wavenumbers for the same parameters and different initial 
conditions). Both systems can demonstrate completely syn¬ 
chronized chaotic state {{a.),D and (b),purple (dark) region) 
complete spatial incoherence {{au),E and (b),blue (light) area). 
Parameters: a = 1.4, p = 0.3, N = 1000. 

By an analytical argument similar to [6| |18| (Appendix) 
one can calculate the critical value of the coupling strength 
CTc (horizontal line in Fig. 0, below which the spatial wave 
profile (fc = 1, period-2 in time) in the coherence tongues 
becomes discontinuous, and a chimera state appears: 

(Tc = 1 - l/\fx{x*,y*) -^1 = 1- l/|2aa:* -h P\. (4) 

where fx is the derivative of / with respect to x, and 
{x*,y*) are fixed points of the map. The analytical re¬ 
sult CTc ~ 0.52 agrees very well with the numerically ob¬ 
tained CTc in Fig.0 and with the snapshots (left column) 
and the space-time plots (right column) shown in Fig. 
for the transition from coherence (Fig. 0) to incoherence 
with decreasing coupling strength. In the tongue C of 
Fig. |3|i for cr < (Jc two spatially incoherent domains ap¬ 
pear (Fig.0) and become wider (Fig. 0-e) with decreas¬ 
ing coupling strength a and, eventually, the dynamics be¬ 
comes completely incoherent (Fig. 0). This scenario rep¬ 
resents a “coherence-incoherence” transition. Similar re¬ 
sults were obtained in 0 for ring networks of logistic maps 
and Rossler systems. 

Figure[0) illustrates the main regimes of nonlocally cou¬ 
pled Lozi maps in the (r, cr) parameter plane. It shows 
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Fig. 4: Coherence-incoherence bifurcation for coupled Henon 
maps for coupling radius r = 0.3 and decreasing coupling 
strength (a) cr = 0.53, (b) cr = 0.47, (c) a — 0.45, (d) cr = 0.32, 
(e) cr = 0.29, (f) a = 0.2). Snapshots at time t = 5000 
(left columns) and space-time plots (right columns) are shown. 
Other parameters are as in Fig. 0. 

a completely synchronized chaotic state (Fig. [3>, purple 
(dark) region), a completely incoherent state (Fig.|0), blue 
(light) region) and regimes of stationary patterns and trav¬ 
elling spatial profiles (Fig. [1>. white region). Similar re¬ 
sults including travelling waves were obtained for the ring 
of coupled Lorenz systems [^. It should be noted that 
chimera states have not been obtained in the ring of Lozi 
maps. The critical value of the coupling strength for cou¬ 
pled Lozi maps, calculated by a similar argument as above 
under the hypothesis of spatial wavenumber fc = 1 and 
period-2 in time (see Appendix), is for cr* < 0 

cTc = 1 - l/|a -/3|. (5) 

which yields CTc « 0.09. However, this value is so small 
that the coupling cannot induce a coherent spatial pro¬ 
file with regular time-period-2 dynamics. Therefore the 
transition to partial spatial coherence does not exist. 

A more detailed analysis of the snapshots and space- 
time scenarios of coupled Lozi maps as protoptype of a 
hyperbolic system reveals completely different scenarios 
from complete coherence to complete incoherence with de¬ 
creasing coupling strength, as shown in Fig. 0 for the case 
of random initial conditions and in Fig.0for specially pre¬ 
pared initial conditions corresponding to a fc = 1 spatial 
wave profile. In Fig. 0r -b the profiles are coherent (com¬ 
pletely synchronized), in Fig. 0; some solitary oscillators 
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desynchronize, whose density increases (Fig. [^), until in 
Fig.[5j3 the system is completely incoherent (desynchro¬ 
nized). In Fig. §r-c the profiles are again coherent, but 
Fig.[ 6 ji,e exhibits discontinuous spatial profiles, however, 
the incoherent parts are not confined to two incoherent do¬ 
mains located between two coherent domains as in nonhy- 
perbolic systems, i.e., no chimera. Rather the whole spa¬ 
tial prohle is incoherent, but the disorder is spread around 
the upper and lower branches. This qualitatively differ¬ 
ent coherence-incoherence bifurcation occurs at a critical 
coupling strength which corresponds approximately to the 
value found analytically in Eq.([^. Finally, in Fig. com¬ 
plete incoherence is reached. 
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Fig. 6: Coherence-incoherence transition for coupled Lozi maps 
for coupling radius r = 0.147 and decreasing coupling strength 
(a) a = 0.92, (b) a = 0.46, (c) cr = 0.18, (d) u = 0.10, 
(e) a = 0.03, (f) a = 0.02. Snapshots at time t = 5000 
(left columns) and space-time plots (right columns) are shown. 
Other parameters are as in Fig. Initial conditions were 
prepared as A: = 1 spatial wave profile. 


system 


Fig. 5: Coherence-incoherence transition for coupled Lozi maps 
for coupling radius r = 0.147 and decreasing coupling strength 
(a) a = 0.94, (b) a = 0.44, (c) a = 0.22, (d) cr = 0.14, (e) 
a = 0.001). Snapshots at time t = 5000 (left columns) and 
space-time plots (right columns) are shown. Other parameters 
are as in Fig. 03. Initial conditions were randomly distributed 
in square xoi G (—0.5; 0.5), yoi £ (—0.5; 0.5). 


Coherence-incoherence transition and chimera 
states in nonlocally conpled Lorenz systems. — 

The absence of chimera states in a ring of coupled Lozi 
maps allows us to advance the following hypothesis. The 
chimera state cannot be obtained in a network of hyper¬ 
bolic systems (both discrete and time-continuous). The 
Lorenz system is hyperbolic for standard parameters. In¬ 
deed, a detailed examination of a ring with nonlocally cou¬ 
pled Lorenz systems 33 confirms this hypothesis. This 


network exhibits complete synchronization, completely in¬ 
coherent states, and stationary patterns or coherent trav¬ 
elling waves but no chimera states. All these regimes were 
also obtained in the network of Lozi maps, but we have 
not found chimera states. 

To test our hypothesis, we have examined the Lorenz 


= —^{x—y) y = —xz+px—y z = —bz+xy, ( 6 ) 


where x, y, z are the dynamic variables, and 7 , b, p are pa¬ 
rameters, in a range of parameters where the Lorenz sys¬ 
tem is nonhyperbolic, i.e., it has a spiral quasiattractor 


37 39 


Fig. [7| il lustrates the bifurcation diagram of the Lorenz 
system ® in the (p, 7 ) parameter plane for fixed parame¬ 
ter 6 = 8/3 [^[^ . The hatched area in the (p, 7 ) param¬ 
eter space shows the hyperbolic regime where the Lorenz 
attractor is found. If 7 = 10 is fixed, then a transition to 
the regime of the nonhyperbolic attractor of the saddle- 
focus type occurs when the bifurcation line Ij, is crossed 
from left to right. 

We simulate a ring of nonlocally coupled Lorenz sys¬ 
tems 40 with parameters p = 220, 7 = 10, b = 8/3, 
corresponding to the nonhyperbolic regime. Indeed, we 
find a coherence-incoherence scenario via chimera states 
composed of one incoherent and one coherent domain as 
shown in Fig.[^ The reason why the chimera does not con¬ 
sist of two incoherent and two coherent domains is that it 
bifurcates from the k = 0 (homogeneous) coherent profile, 
and not from the k = 1 coherent profile as in Fig.4. Such 
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Fig. 7: Bifurcation diagram of the Lorenz system Eq. § show¬ 
ing the hyperbolic regime of the Lorenz attractor (hatched) 
and the nonhyperbolic attractor for h = 8/3. 


different types of chimeras have also been noted in the Van 
der Pol system [I^. 



of such complex systems. 
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Appendix: Critical coupling strength. — In 

the following, we analyze the bifurcation scenario from 
complete coherence to complete incoherence via chimera 
states, and derive an analytic estimate of the critical cou¬ 
pling strength for the onset of chimera states for the 
Henon and the Lozi map, respectively. 

From a geometrical point of view and in the thermody¬ 
namic limit A —>■ oo, where r = P/N is fixed, coherent 
solutions x* approach a smooth profile x*{^), where ^ is 
the spatial variable, and the spatially continuous version 
of Eq. ^ is given by 

A transition from coherence to partial incoherence 
(chimera state) occurs if the respective profile x*{^),y*{^) 
becomes discontinuous at some points ^ of the ring 5^. 

Let us consider a solution of system ([^ with wave num¬ 
ber k = 1 (i.e., wavelength equals system size) and period- 
2 dynamics in time. Hence we can reduce the dynamics by 
even and odd time steps and x^{^),y^{^), re¬ 

spectively. This leads to (^) = (1 —cr)/(A (^), j/^ (^)) + 
^ I^-r f i^^(d),yHv))dv, y^~HO = Idx^iO with j = 
0,1. Taking the spatial derivative yields 


Fig. 8 : Chimera state in a ring of A = 300 nonlocally coupled 
Lorenz systems in the nonhyperbolic regime. Snapshots at time 
t = 5000 (left column) and space-time plots (right column) are 
shown. Parameters: b = 8/3, 7 = 10, p = 220, r = 0.22. 
Decreasing from top to bottom, cr = 0.4, cr = 0.35 and o = 
0 . 1 . (a) - complete synchronization, (b) - chimera state, (c) - 
spatial incoherence. 


Conclusions. — We have hypothesized that chimera 
states can be obtained in networks of chaotic nonhyper¬ 
bolic systems and that they cannot be found in the net¬ 
works of robust hyperbolic systems. This hypothesis has 
been confirmed by numerical simulations of a ring of non¬ 
locally coupled time-discrete and time-continuous systems 
both for hyperbolic and nonhyperbolic cases. 

Our results may be applied to a variety of diverse sys¬ 
tems. On one hand, chaotic behaviour can be found in nu¬ 
merous natural systems such as weather and climate, and 
chaotic dynamics plays an important role in many fields of 
science and engineering ranging from physics, chemistry, 
biology to economics and sociology. On the other hand, 
there has been growing interest in chimera patterns with a 
wide scope of applications. Therefore, disclosing the uni¬ 
versal mechanism of the emergence of chimeras in chaotic 
systems is a vital step in understanding the functionality 


2^5 ■^(0 = (l-c^) fxi{.x^ ,y^)x\ +fyj{a 


,y^)y\ 


+ ^ [/ {x^+ x),y^+ r)) - f {x^{,i-r),y^{i-r))] 

yl{0=H-\0- ( 8 ) 


where fx{x,y) = dxf(x,y) etc. At the point ^ where the 
smooth profile breaks up, the spatial derivative becomes 
infinite. Considering that |a;||,|y|| diverge to infinity, we 
can neglect the coupling term on the right-hand side of 
Eq. ([^. The main contribution comes from the first term. 
Using fy{x,y) = 1 we obtain: 



fxd 

^,y^)x^^ +/3x^ ^ 

(9) 

Multiplying the equations 

for even and 

odd 

time steps (j = 0 , 1 ) we 

obtain x'^x^ = 

(1 - 

0 -)^ fx^{x^,y^)x\+ Px^ 

fxo{x 

°,y°)x° + /3x^ . 

Now 


assume the symmetry of the spatial profiles at the 
turning points (where the slope m becomes infinite) 
a;^ = —= m, which yields the following condition 


1 = (1 - af [fx^{.x^,y^) - /3] [fxo{.x°,y°) - j3] . (10) 


In order to obtain an analytic approximation for the crit¬ 
ical coupling strength CTc at the onset of the chimera 
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states, we assume that the points ^ where the smooth 
profile breaks up and the spatial derivative becomes in¬ 
finite are hxed points {x*,y*) of the map. This yields 


= X*, = y* and Eq.(|10[) becomes 


1 = (1 — af [fx{x*,y*) — j3] . This leads to 


(Tp = 1 ± 


\fx{x*,y*) - /3| ■ 


( 11 ) 


where the minus sign should be chosen since the lower 
value of CTc represents the threshold where the smooth 
profile breaks up with decreasing coupling strength cr, 
and the spatial coherence is lost. 


For the Henon map f{x, y) = 1 — ax^ + y the derivative 
is equal to fx = —2q;x*, and we obtain CTc = 1 — \ 2ax ‘‘+ i3 \ ■ 
The fixed points of the Henon map for a = 1.4, (3 = 0.3 
are approximately x* = 0.63 and x* = —1.13. Using x* = 
0.63 we find CTc ~ 0.52. 

For the Lozi map f{x,y) = 1 — a\x\ + y the deriva¬ 
tive is equal to fx = —a for x* > 0 , and fx = a for 
X* < 0, and we obtain for x* < 0: (Jc = 1 — 

The hxed points of the Lozi map for a = 1.4, (3 = 0.3 
are approximately x* = 0.48 and x* = —1.43. Using 
X* < 0 we hnd CTc ~ 0.09, which is so small that no 
period -2 dynamics exists, since the dynamics is chaotic 
(similar to the uncoupled maps). Hence the above 
argument is not valid, and a transition to partial in¬ 
coherence (chimera state) of the above type does not exist. 


From the numerical simulations for the Henon map we 
obtain that the transition from coherence to incoherence 
in system ^ occurs for a coupling strength a close to 0.5. 
The deviation from our approximation is due to the hnite 
number of nodes in the simulations, while the analytics 
assumes the continuum limit. 
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